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Csl ; Abstract 

QJ ' Mesoscopic theory for soft-matter systems that combines density functional and statistical field 



X 



theory is derived from the microscopic theory by a systematic coarse-graining procedure. Grand- 
thermodynamic potential functional for hard spherical particles that interact with (solvent medi- 



ttn ' ated) spherically-symmetric potentials of arbitrary form is a sum of two terms. In the first term 

O ' 

■ microscopic length-scale fluctuations are included, and the second term is the contribution associ- 



ated with mesoscopic length-scale fluctuations. In the approximate theory the flrst term has the 
form of the density functional in the local density approximation, whereas the second term has the 



Q ' form known from the field theory and depends on the pair correlation function for which a pair of 

equations similar to the Ornstein-Zernicke equation with a particular closure is obtained. For weak 
1^ , ordering the theory can be reduced to the Brazovskii field theory with the effective Hamiltonian 

I having the form of the grand-potential functional in the local density approximation. 

On ■ 

^SJ ' Within the framework of this theory we obtain and discuss the A-line and the universal sequence 



of phases: disordered, bcc, hexagonal, lamellar, inverted hexagonal, inverted bcc, disordered, for 



o 

00 

I increasing density well below the close-packing density. The sequence of phases agrees with ex- 



perimental observations and with simulations of many self-assembling systems. In addition to the 
above phases, more complex phases may appear depending on the interaction potentials. For a 
particular form of the short-range attraction long-range repulsion potential we find the bicontinu- 
ous gyroid phase (Ia3d symmetry) that may be related to a network forming cluster of colloids in 
a mixture of colloids and nonadsorbing polymers. 



I. INTRODUCTION 



Various soft-matter systems, including surfactant solutions, globular proteins, colloids in 
different solvents, colloid-polymer or star polymer - linear polymer mixtures, exhibit self- 
assembly into various structures. In particular, clusters or micelles of various shape and size, 
including branched networks can be formed. The clusters or micelles may exhibit ordering 
into different periodic phases for a range of volume fractions, including quite dilute systems 

Phase-separation competes in such systems with formation of 
lyotropic liquid crystalline phases. Transitions between soft crystalline phases having differ- 
ent symmetries, phenomena such as re-entrant meltingjlO, 11, 12, 13|, and finally formation 
of hard crystals are observed when the volume fraction increases. Theoretical studies of 
such systems can be based on liquid theories such as SCOZA, Landau-type theories or den- 
sity functional theories (DFT). Unfortunately, the SCOZA cannot describe the formation of 
inhomogeneous microphases 14l |. 

Landau-Ginzburg- Wilson (LGW) Q theory turned out to be very successful in de- 
scribing universal properties of critical phenomena associated with transitions between uni- 
form phases. Landau-Brazovskii (LB) 13, Q, [3] theory on the other hand predicts universal 
properties of transitions between uniform and periodically ordered lyotropic liquid crystals. 
The statistical field-theoretic methods are very powerful in determining effects of long-range 
correlations between fluctuations that lead to ordering - either to phase separation or to 
formation of various ordered structures on the nanometer length scale. The effective Hamil- 
tonians in Landau theories are given in terms of phenomenological parameters whose pre- 
cise relation with thermodynamic and material properties is irrelevant for the universal, i.e. 
substance independent properties. Unfortunately, the general, abstract theory has limited 
predicting power for particular systems. 



On the other hand, in the very successful density functional theory 20|] the short-range 
structure is taken into account quite precisely, whereas long- and intermediate range scale 
fluctuations are taken into account in a very crude mean-field approximation in most appli- 
cations of the DFT. The DFT and LGW or LB theories are complementary in the treatment 
of the short- and long range correlations. It is thus desirable to develop DFT with the form 
of the grand thermodynamic potential functional that includes the contribution associated 
with the mesoscopic scale fluctuations. In order to develop such a theory one should per- 
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form a systematic coarse-graining procedure with controllable accuracy, and derive effective 
Hamiltonians from the microscopic ones. Nonuniversal properties are correctly described 
within the collective variables (CV) 
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23l ] and hierarchical reference (HRT) 
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theories. However, the HRT is restricted to critical phenomena, and the formal structure of 
these theories is rather complex. 

Here we propose an alternative approach, resulting in a density functional theory with a 
rather simple structure. The theory allows for including in the grand potential the contribu- 
tion associated with long-range correlations between fluctuations. The latter contribution 
can be calculated within field-theoretic methods. Within the framework of our theory it is 
possible to determine phase transitions between different soft-crystalline and uniform phases. 
In special cases of phase separation or weak ordering into soft crystals the theory reduces 
to the standard LGW or LB theories respectively, with the coupling constants expressed in 
terms of density, temperature and the interaction potential. Within the present approach 
phase diagrams in terms of density and temperature, rather than abstract phenomenological 
parameters can be obtained, and validity of the LGW or LB theories in particular systems for 
given thermodynamic conditions can be verified. This kind of approach was applied already 



to highly charged colloids 
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261 and to ionic systems 
The derivation of the theory is described in sec. II. In sec. Ha the mesoscopic state 
(mesostate) is defined, and in sec. lib the probability distribution for the mesostates is 
derived from the microscopic theory. In sec. Ilia the correlation functions for the meso- 
scopic densities are introduced and their relation with the microscopic correlation functions 
is discussed. Vertex functions and the grand-potential density functional are introduced 
in sec. Illb. Self-consistent equations for the two-point vertex functions (related to direct 
correlation functions) are derived in sec. IIIc. Periodic structures are considered in sec. 
Illd, where the approximate expression for the grand potential is also given. In sec. IV the 
general framework of the theory derived in the preceding sections is applied to a particular 
approximation for the probability distribution for the mesostates, related to the local den- 
sity approximation. The self-consistent equations for the two-point functions reduce in this 
approximation to simpler forms given in sec. IV a, and the approximation for the grand 
potential is given in sec. IVb. The relation between the present theory and the Landau-type 
theories is discussed in sec.IVd. Explicit results for the case of weak ordering are briefly 
described in sec. V. In sec. Va we study the A-line in the Brazovskii theory. In sec. Vb we 



3 



limit ourselves to universal features of the phase diagrams obtained in the simplest one-shell 
mean field (MF) approximation, and show the universal sequence of phases: disordered, bcc, 
hexagonal, lamellar, inverted hexagonal, inverted bcc and disordered for increasing density 
of particles (well below the close-packing density). The details of the phase diagrams can 
be obtained beyond the one-shell MF approximation. However, these details depend on the 
shape of the interaction potential. Studies of particular systems go beyond the scope of this 



work and will be described elsewhere 



361 ]. For an illustration we consider a particular form 



of so called short-range attraction long-ran 



considerable interest recently 
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sion (SA 
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R) potential that has drawn 
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4J|. In sec. Vc we quote 



results for a particular form of the SALR potential in the two-shell MF approximation, 
and show stability of the bicontinuous gyroid phase between the hexagonal and the lamel- 
lar phases. Thisphase may be related to the network-forming cluster of colloids observed 
experimentally [7|. 



II. COARSE GRAINING 



Let us consider hard spherical objects that interact with arbitrary spherically-symmetric 
potentials. The approach can be generalized to hard object of different shapes, but it is easier 
to fix attention on spheres. A microstate is given by the sequence {ra}a=i...N describing the 
positions of the centers of spheres. In addition to the microscopic density, 

/)(r,{rj):=^<5(r-r„) (1) 

a 

we consider the microscopic volume fraction 

f,{r,{r^}):=Y^eia-\r-r^\) (2) 

a 

where a is the diameter of the hard sphere. The energy in the microstate {'ra}a=i...N is given 
by 

E[M] = T.nK-r,\), (3) 

a>l3 

where l^(|r — r'|) is the pair interaction potential. 
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A. Mesoscopic density and mesoscopic volume fraction 



Let us choose the mesoscopic length scale R> a /2 and consider spheres Sr{v) of radius 
R and centers at r that cover the whole volume V of the system (for the bulk system we 
assume periodic boundary conditions). We define the mesoscopic density and the mesoscopic 
volume fraction at r by 



respectively, where Vs = AttR^/S. For brevity we shall use the notation J,. = J dr, indicating 
the integration region 5* by /r/gg when necessary. 

For given mesoscopic length scale R, and given microstate {r^}, ?7(r) is a continuous field 
such that r]{r) < rjcp for all r, where rjcp is the close-packing volume fraction. Moreover, 
the gradient of r] is small, IVt^I < 1/R . The mesoscopic density has discontinuities. The 
steps are n/Vs, where n is the difference between the number of the centers of hard spheres 
included in Sr{t) and Sr^t + dr). 

The microstate {fq} for which (jll) (or ([5])) holds for all r is called compatible with the 
field p(r) (or //(r)) for given R. In order to indicate that the microstate {tq,} is compatible 
with the field p(r) (or ^^(r)), the notation {r^} G {p{r),R} (or {r^} G {ri{r),R}) will be 
used. The set of all microstates compatible for given mesoscopic length scale R with the 
field p(r) (or ^^(r)) is called a mesostate at the length scale R. 

For given mesoscopic length scale R the set of all microstates can be split into disjoint 
subsets, such that each subset contains all microstates compatible with a particular field 
p(r) (or 77(r)) according to the definition (jH (or (JS])), and no other microstates. 

The mesostate gives less detailed information about the state of the system than the 
microstate, but more detailed than the macrostate. The mesostates depend on the chosen 
length scale (see Fig.l). For R ^ oo the mesostates become identical with the macrostates. 
This is because so defined mesostate is characterized by the number density (or volume 
fraction of particles) in the whole system. (Periodic boundary conditions are assumed). 
On the other hand, when R = a/ 2, Eq.(l5]) defines the one-to-one relation between the 
microstates and the mesostates. To prove the above let us consider two different microstates 
{fq,} and {r'^}. If they are different, then a center of at least one sphere in the first microstate 
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FIG. 1: The mesoscopic volume-fraction rj defined in Eq.([5]) (dashed hnes) and pvQ (sohd hnes), 
where p is defined in Eq.Q and vq = ttit^/G, for the microstate in which a single hard sphere of a 
radius a/2 is located at r = 0. Top panel: 2R/a = 1. Central panel: 2R/(t = 2. Bottom panel: 
2R/a = 5. In this particularly simple case the fields ([5]) and (H]) are functions of the distance r 
from the center of the hard sphere. In each case /j. 7y(r) = vq. The distance r is in a/2 units, r] and 
pvQ are dimensionless. 
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does not coincide with a center of any sphere in the second microstate. Let us assume that 
ri 7^ y' 13 for all (3 = 1,...,N. Let us calculate ([5]) for r = ri. For the first microstate we 
get ?7(r) = 1, for the second microstate ?7(r) ^ 1. Thus, different microstates cannot be 
compatible with the same field defined in 

The choice of R depends on the problem under consideration, in particular on the length 
scale A characteristic for ordering, and is determined by the Hamiltonian. In order to 
describe the ordering on the length scale A, one should choose a/2 < R < A/2. 

B. Probability distribution for the mesostates 

Let us calculate the probability density that the mesostate {p(r), R} (or {//(r), R}) occurs 
in the system spontaneously. This is equal to the probability that any microstate compatible 
with p(r) (or ^^(r)) occurs. We derive the expressions in terms of p; in terms of rj the theory 
has the same formal structure. In an open system in contact with the thermostat the 
probability density of the mesostate {p, R} is given by 

p[p] = l[ e-^(^-'^X''W) (6) 

where H is the microscopic Hamiltonian, /5 = 1/(A;bT) and /{r„}g{p,_R} is the symbolic nota- 
tion for the integration over all microstates compatible with p. p and T are the chemical 
potential and temperature respectively, and (3 = l/ksT, with ks denoting the Boltzmann 
constant. Finally, 

J{rc.} J J{r^}e{p,R} 

The functional integral / Dp in ([7]) is over all mesostates {p,R}, which is indicated by the 
prime. 

Fixing the mesostate in the system is equivalent to the constraint on the microstates of 
the form (jlj). In the presence of the constraint {p,R} the grand potential is denoted by 
^co[p] and is given by 

i{r„}e{p,i?} 

Thus, the probability density of a spontaneous occurrence of the mesostate {p, R} (equal to 
the probability density that any of the microstates compatible with p occurs) is given by 

P[P] = (9) 
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where 




(10) 



We obtain a mesoscopic theory with the same structure as the standard statistical me- 
chanics. The integration over all microstates is replaced in ( llOp by the integration over 
all mesostates. The Hamiltonian is replaced in ([9]) and ( ITOl) by the grand potential in the 
presence of the constraint of compatibility with the given mesostate that is imposed on 
the microstates. The above formulas are exact. So far we just rearranged the summation 
over microstates. In ffTOl) the integration over microstates compatible with each mesostate 
is included in Qco, and then we perform a summation over all mesostates. The reason for 
doing so is the possibility of treating the summation over the mesostates and over the mi- 
crostates compatible with a particular mesostate on different levels of approximation. In a 
similar way one can define a mesoscopic theory for the mesoscopic volume fraction. Both 
approaches are equivalent. The advantage of the density field is the fact that the microscopic 
and functional-density theories are based on the density rather than on the volume fraction, 
and the form of Qco should be determined within microscopic theories. The advantage of 
the volume fraction is the fact that it is a continuous field. 

Note that in the mesoscopic theory the mesostate {po, R} that corresponds to the global 
minimum of fico[p] is analogous to the ground state in the microscopic theory (similar role 
plays riQ corresponding to the minimum of f2co[r7]). This is because the ground state is the 
microstate that corresponds to the global minimum of the Hamiltonian. In this context 
the important property of the ground state is the fact that it corresponds to the most 
probable single microstate. Likewise {pq, R} is the single mesostate that occurs with the 
highest probability, because iQ assumes a maximum when fico[p] assumes a minimum. 

III. GRAND POTENTIAL AND THE CORRELATION FUNCTIONS FOR THE 
MESOSCOPIC DENSITIES 

The grand potential in the system subject to the constraint for the mesoscopic density 
distribution p(r) can be written in the form 



'CO 



U-TS- pN, 



(11) 
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where U, S, N are the internal energy, entropy and the number of molecules respectively in 
the system with the constraint (jl]) or ([5]) imposed on the microscopic densities or the volume 
fractions respectively. U is given by the well known expression 

U[P*] = \l I ^-(ri - r2)p*(ri)p*(r2), (12) 

where 

Vcoi'Oi - Y2) = V{ri2)gco{ri - r2), (13) 

^12 = ki — r2|, and gco{^i~T^2) is the microscopic pair correlation function in the system with 
the constraint (jll) imposed on the microscopic states. From now on we consider dimensionless 
density 

p* = pa\ (14) 

In the case of the considered systems (no internal degrees of freedom of the particles) 
the entropy S satisfies the relation —TS = Fh, where is the free-energy of the reference 
hard-sphere system with the constraint (j4]) or (j5]) imposed on the microscopic densities or 
the microscopic volume fractions respectively. 

A. Correlation functions for the mesoscopic densities and their generating func- 
tional 

Let us introduce an external field J(r) and the grand-thermodynamic potential functional 
Q[f3J] := -ksTlog [ J Dp*e-^[^-[''*l-X^W^*Wl]. (15) 

The generating functional for the (connected) correlation functions for the mesoscopic 
densities is — /5fi[/5J], and 

^'^^-P ^'-^^ -5(/5J(rO)...5(/3J(r„))- ^'^^ 
In particular, the above gives 

where (jlj), ([8]) and ([7]) were used. We use the same notation (...) for the microscopic and 
the mesoscopic average. The average value of the mesoscopic density at the point r is the 



average microscopic density integrated over the sphere Sr{y) of the radius i?, and divided 
by its volume Vs- Clearly, (p*(r)) depends on R. If, however, (p*(r)) is independent of R for 
some range of cr/2 < R < A/2, then this function gives information about actual ordering 
on the length scale A, because the average value of the mesoscopic density is the microscopic 
density averaged over regions smaller than the length scale on which the ordering occurs. 
For the correlation function for the mesoscopic density we introduce the notation 

g2{r,-v,) = {p\vr)p\v,)Y°\ (18) 

where here and below the superscript con denotes 

{A{p*{v,))B{p*{v,))r- = {A{p*{r,))B{p*{r,))) - {A{p*{r,))){B{p*{r,))), (19) 

and where 

l^y,, /Z^pe-/^["-'^1-X^(-)^*(-)]p-(r,)p*(r,) 



(p(r')p(r")). 

The correlation function for the mesoscopic density at ri and r2 is the correlation function 
for the microscopic density at the points r' G SfSji) and r" G Sb^Ji)-, integrated over the 
spheres Sr{yi) and Sr{y2) respectively. While the microscopic pair distribution function 
g^ir' - r") = ^(r' - r")p(rXr") = (p(r')p(r")) - (p(r'))5(r' - r") vanishes for |r' - r"| < (X, 
the mesoscopic quantity defined by Eqs. fll8l) - fl20l) does not vanish for ri = Y2 when R > cr/2 
(see Fig. 2). The larger R is, the larger is (p(ri)p(r2)) — (p(ri))5(ri — r2) for ri — > r2. 

B. Grand-potential functional of the mesoscopic density, and the vertex functions 

Let us introduce the Legendre transform 

(3F[P*] ■■= m/3J] + [pj{v)p*{v) (21) 

where 

fir) - ^(-^^^ (22) 
is the dimensionless average field for given J(r). The equation of state takes the form 
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FIG. 2: Microscopic configuration contributing to the correlation function ^2(ri — for ri = r2 
in the case of R > a/2. The black circles represent particles, and the open circle represents the 
coinciding spheres Sji{ri) = <S'_R(r2), over which the microscopic correlations are averaged. Note 
that the black spheres do not overlap, and their centers are both included in Sji. Thus, the 
microscopic correlation function for this microstate does not vanish. 

In general p* may differ from any mesostate defined by (jlj). Let the functional Qco, 
defined for the mesostates, be extended beyond the set of the mesostates. Let the extension 
be defined in (11 II) on the Hilbert space of fields / that fulfill the restrictions following 
from the properties of the mesostates, and let us keep the notation Qco for this extension. 
The key restriction on the mesoscopic volume fraction is the gradient, whose upper limit 
is a'^/{2RY < 1/{2R). The mesoscopic number density has discontinuities, with the steps 
~ / B? decreasing for increasing R/a (see Fig.l), and the separation between the steps is 
~ o". In the approximate theory we shall consider the Hilbert space of fields such that in 
Fourier representation /(k) = for the wavenumbers k > n/R. Here and below /(k) is the 
Fourier transform of /(r). The mesoscopic volume fraction and density are bounded from 
above, and close-packing is the natural limit for their values. However, the fields p* (k) with 
large values yield large values of Qco (see f lTT]) and f|T2|) ). and in turn small values of the 
Boltzmann factor (Q, therefore inclusion of such fields in the Hilbert space should not have 
large effect on the results. Selfconsistency of the approach requires that the average density 
obtained within this theory is bounded from above by the close-packing density. 

We can define 



where p* is given in (!22l) and C^°(ri, ...r„|p*] is the n—th functional derivative of fico[p*] at 
p* = p*. The second equality holds for the extended functional, for which the derivative can 




(24) 
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be defined. By definition 



0. Then from ffTSD and fl2Tl) we obtain 



l3n[l3J] = -I3n,,[p*] + / l3J{v)p*{v) + log / Dct^e 



(25) 



and 



PF[p*] = pn,,[p*] - log / D0e 



(26) 



Eq. f l26p defines the functional in which the microscopic scale fluctuations, with frozen 
fluctuations of the mesoscopic density over the length scales larger than R, are included in the 
first term, and the mesoscopic scale fluctuations are included in the last term. The expression 
( !26|) would be exact, if the functional integration were restricted to the mesostates. However, 
the powerful methods of the functional analysis could not be applied. There are no easy 
ways of estimating the error associated with inclusion in the fluctuation contribution fields 
that do not represent the actual mesostates. The Eq.( l26i) should serve as a starting point 
for various approximate theories. The idea is to apply approximate microscopic theories to 
the first term, and field-theoretic approaches to the second term. 



C. Equations for the correlation functions 



Note that from fl26l) it follows that the vertex functions (related to the direct correlation 
functions) defined by 



(27) 



5p*(ri)...5p*(r„) 

consist of two terms: the first one is the contribution from the fluctuations on the microscopic 
length scale (< R) with frozen fluctuations on the mesoscopic length scale, the second one 
is the contribution from the fluctuations on the mesoscopic length scale {> R). 
The average density for J = satisfies the equation (see Eqs. (|25il and (12^ ) 



0, 



(5p*(r) 6p*{v) 

where the averaging is over the fields with the probability oc exp(— /3ifjiuc[p* 
two-point function from fl27|) and fl26l) we obtain the equation 



2C2(ri,r2)=Cr(ri,r2) + 



6p*{Ti)6p*{Y2) 



)-(- 



5H 



flue 



^6p*{ri] 



[r'))C2(r',r2) + 



6p*{ri) 6p*{T2 



'(5p*(r2) 



r'))C2(r',rO 



(28) 
For the 

(29) 
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In calculating the functional derivative of F (Eg. (126!) ) we used the equality 5J {vi) / 5p* {r2) = 
C2(ri,r2) (see ((23])). From ^ we have 



5p*(r)...5p*(r(™)) ^Jr^ Jr„ n! 

Higher order vertex functions can be obtained from Eqs. fl27|) and fl26l) in a similar way. As 
a result, a hierarchy of equations where the vertex functions are expressed in terms of the 
many-body correlation functions is obtained. 

Eqs. ( |28l) - (!30l) allow for including the fluctuation contribution to the vertex functions on 
different levels of approximation. In the mean-field approximation the fluctuation contribu- 
tion is just neglected. In the lowest-order nontrivial approximation only terms proportional 
to the pair-correlation function (0(r)0(r')) are included in fl28l) - (!29l) . by truncating the ex- 
pansion in fl5Ul) at = 2. Then Eqs. (1281) and fl2^ assume the approximate forms 

6p*{r) 

and 



+ / / ^2(ri,r2)C3^°(ri,r2,r) = (31) 



2C2(ri,r2) = 3Cr(ri, r^) + ^ G^ir, r') f^"^""' ""'^ ""^^ - C-(r, rO^" (r^ r^)] . (32) 
From fl23l) and ([T6l) we obtain the analog of the Ornstein-Zernicke equation 



C2(ri,r2) 



6p*{r2)6p*{ri) 
or in the equivalent form 



Sp*ir2) 
W(ri), 



-1 



e^2(ri,r2)-\ (33) 



/ C2(ri,r2)e?2(r2,r3) =<5(ri-r3). (34) 

We assume that the form of flco is known from the microscopic theory, and consider it as an 
input to the mesoscopic theory. Selfconsistent solutions of Eqs. fl32|) - fl34|) yield the two-point 
vertex and correlation functions, C2 and Q2 respectively, for given forms of Note that 
Eq. (l32l) plays a role of the closure to the OZ equation fl34l) . 

Eq. (pTj) is the minimum condition for the grand potential. Since there may exist several 
local minima, the solution corresponds to a stable or to a metastable phase when the grand 
potential assumes the global or the local minimum. The solution of Eq.( l3Tl) corresponding 
to the global minimum gives the average density for given /z and T in the lowest nontrivial 
order beyond MF. 
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D. Periodic structures; general case 

Let us consider periodic density profiles 

p*ir)=pl + ^r) (35) 

where 

$(r + P) = $(r) (36) 

and P = J2i ^iPi where Pi are the vectors connecting the centers of the nearest-neighbor 
unit cells and rii are integer. The pg is the space-averaged density, i.e. 

m = 0, (37) 

reVu 

where Vu is the unit cell of the periodic structure, whose volume is denoted by Vu- In the 
case of periodic structures 

C2(ri + P, r2 + P) = C2(ri, rs) = C2(Ar|r2) (38) 

where Ar = ri — r2 G i?'^ and r2 G V^. 
Let us consider the Gaussian functional 

^g[p*, 0] = ^ / / 0(ri)C2(ri - r2|r2)0(r2) = ^ / / 0(k)C2(k, k + k')0(k'). (39) 

Here and below we use the simplified notation = / c/k/(27r)'^. For periodic structures we 
have (see Appendix A) 

C2(k, k + k') = 5(k + k')C'2(k) (40) 

where 

(72(k) = / C2(Ar)e^^■^^ (41) 



and 



C2(Ar) = / C2(Ar|r2) (42) 



is the vertex function averaged over the unit cell. For periodic structures the functional ([39 
can be rewritten in the equivalent form 



^g[p*,0] = \ /^0(k)C2(k)0(-k). (43) 
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From the analog of the Ornstein-Zernicke equation (133|1 we obtain the relation 

C'2(k)G'2(k) = 1 (44) 

between the vertex function averaged over the unit cell, and 

G2(k)=/ / e^'^-^'-^-^^^e^slri.ra) = / e'''<''--''-^G^{v^ - v^). (45) 

The function 6*2 (^i — r2) describes the correlations between the density fluctuations at the 
distance ri — r2, with the first point position averaged over the unit cell. 

In order to calculate the second term in ( !26l) we follow the standard field-theoretic pro- 
cedure and write 

HfiuAP*, 0] = Hg[p*, 0] + ^W*, 0] (46) 

where Ti^lp*, 0] is the Gaussian functional ( l43l) . Next we make an assumption that A Ti [p * , ij)] 
can be treated as a small perturbation. When such an assumption is valid, we obtain 19l.l31| 

(3n[-p*] ^ (3n,,[p*] - log I D0e-^^« + {PAn)G + 0{{pAn)l). (47) 

where {■■■)g denotes averaging with the Gaussian Boltzmann factor e~^'^^ . Note that the 
Gaussian functional integrals can be calculated analytically. Eqs. ( |32l) - (H71) allow for 
calculation of the fluctuation contribution to the grand potential in the lowest nontrivial 
order when the form of VL^o is known, and the form of C2 ( fj4Tl) and (H2|) ) is determined 
by a self-consistent solution of Eqs. (l32|) and (!33l) . In general, each contribution to Eq. p7|) 
depends on the mesoscopic length scale i?, but the i?-dependent contributions must cancel 
against each other to yield i?-independent Vt. By minimizing the density functional fH7|) we 
find the equilibrium structure. 



IV. APPROXIMATE THEORY: LOCAL DENSITY APPROXIMATION FOR O^o 

In our theory different levels of approximation for fico are possible. For an illustration 
we demonstrate the mesoscopic theory with the microscopic degrees of freedom considered 
in a very crude approximation. In this approximation calculations are greatly simplified 
and analytical methods can be used. Next we show how the functional can be further 
reduced to the Landau-type form, and discuss in what cases such a reduction is justified. 
The expressions for the phenomenological parameters in the Landau-type theory are given 
in terms of the thermodynamic variables and Vco- 

15 



A. Derivation of the functional Qco in the local density approximation 



In order to derive an approximate form of flco (Eg. (11 II) ) we need an approximate form 
of the microscopic pair correlation function. The microscopic correlation function has the 
limiting behavior 

gco{ri - ra) = for |ri - rs] < cr (48) 
5fco(ri - r2) ^ 1 for |ri - ^ cx). 

If the ordering in the system occurs on the length scale larger than R ^ cr/2, then the 
precise form of gco is not crucial and in the simplest approximation we assume 

(7co(ri-r2) = 0(|ri-r2|-a). (49) 

In order to develop an approximation for the entropy S in the presence of the constraint 
imposed on the density profile, let us consider the number of the microstates associated 
with different positions of the centers of particles included in the sphere Snir) for given 
p(r). We assume that the corresponding contribution to the entropy only weakly depends 
on the mesoscopic density at r' when |r — r'| ^ i?. When this assumption is satisfied, then 
the local density approximation can be applied, and 

F,[p*] = Jj,ip*{r)), (50) 

where fh{p*) is the free-energy density of the hard-sphere system of density p* = pa^. For 
the latter we may assume the Percus-Yevick or the Carnahan-Starling approximation. We 
choose the former case (compressibility route) and assume 

"3^(2 - v) 



(3f,{p*) = p*\n{p*)-p*+p* 



— ln(l — 7]) 



2(1 -r/)2 

In this simple approximation we obtain the i?-independent functional 



(51) 



/3ajp*] = \ f f /?K„(ri2)p*(ri)p*(r2) - / PMp*{t)) - [ /?^p*(r), (52) 

A J n Jr2 Jr 

where ^0(^^12) is defined in Eq.( |T3i) . and for gco{j'i2) we make the assumption ( l49l) . Note 
that in this approximation the dominant contribution to the fluctuation part in (H7j) should 
be independent of R, otherwise the theory is not valid. This is because VL cannot depend on 
the arbitrary length scale R. 

16 



For Fh[p*] given by (l50i) . i.e. for the local density approximation, the functional deriva- 
tives of the functional fico take the forms 



C(ri,...rn|p*] 



5(ri-r2)...(5(r„_i-r„)/5/i")(p*(ri)) for n>3 

(P*(ri))5(ri-r2)+K.o(ri2)) for n = 2 (53) 

/?(tf)(p*(ri))-/i + /r,Ko(ri2)p*(r2)) for n = 1. 

where fjj^\p*{r)) denotes the n-th derivative of fh with respect to its argument, calculated 
at p*(r) = p*(r). In the local density approximation Eq. fl2^ assumes the simpler form. 



(3Hfi^,[p*, 0] = f3Q,,[p* + 0] - f3Q,,[p*] = (54) 
U [ 0(ri)C2^°(ri,r2)0(r2) + / Cr(r)0(r) + J] / ^"^^ (r)-. 

The dominant contribution to the second term in fl26l) comes from small fields 0; moreover, 
fields with large values do not represent the actual mesostates. For fields with small values 
the expansion in (154|) can be truncated. For stability reasons the 0^ term must be included. 
When the expansion in flM|) is truncated at the fourth order term, we obtain Hfi^c of the 
form similar to the effective Hamiltonian in the Landau-type y?^ theories, except that in 
nonuniform phases the coefficients that multiply 0" depend on space position in a way 
determined by the form of fh{p*{r)). 

B. Approximate equations for the correlation functions 

In the local density approximation Eqs. (l3T]) and (!32l) assume the simpler forms 

M,/f(m,,,,).o (55) 

and 

2C2(ri,r2) = 3C2^°(ri,r2) (56) 



-5(ri-r2)6;2(ri,ri)- / / Cnr,,r)g2{r,r')C^"{r' ,r,) 



In this approximation Eqs. ( ]56l) and (IMI) should be solved selfconsistently. Eq. ( l55l) is the 
extremum condition for the grand potential. 
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C. Periodic structures in the local density approximation - case of weak ordering 



Let us consider the functional (15^ for the periodic density profiles (135|1 . and let the 
expansion in be truncated at the fourth order term. We restrict our attention to weak 
ordering, $ ^ Pq. For small $(r) the Taylor expansion 

ft\p*ir)) = ft\Po) + E ' J '^ ^i^r (57) 

m=l 

can be truncated. In the consistent approximation we truncate the above expansion at the 
fourth order term in $. In this approximation Eg. (1541) assumes for periodic density profiles 
(!35|) the form 

f3Hf,^4<l>, 0] = - / 0(k)Cr(A:)0(-k) + / Y: ' V (rr. (58) 

where the prime in the above sum means that n 7^ 2, the explicit expressions for are 
given in Appendix B, and the function C2°{k), defined as in Eq. pTl) . has the explicit form 
(see (dmD) 



VUk) + ff\-pi) + ^^<l> 



(59) 



with 



$2 = 1/ $(r)2. (60) 



From ([56]), (jH]), and (HH) we obtain the equation for C2ik) 



2C,{k) = 3C-{k) + ^Jl^g - 3!i^, (61) 

2 C2[k) 

where we introduced the notation 

Recall that by construction of the mesoscopic theory on the length scale i?, the cutoff ~ vr/i? 
is present in the above integral. Recall also that ^\^G2{k) = G2{0) is not the microscopic 
correlation function at zero distance, but the correlation function for the microscopic density 
at the points r' G 5'ij(r) and r" G Sr{jc), integrated over v' G Sr{y) and r" G Sr{y), as 
discussed in sec.IIIA. The average density for given p and T is the solution of Eq. (l55|l . If 
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there are several solutions, the one corresponding to the global minimum should be chosen. 
In practice it is much easier to choose the average density pg as the independent parameter. 

In the case of weak ordering, i.e. for p*(r) given in Eqs. fl35p and (I57|) with $ ^ pg, the 
truncated Taylor expansion of fh{p*{r)) about p*(r) = pg can be inserted in Eq. flS^ (see 
( |Fr|) for n = 0), and PQco[p*] can be approximated by 

pn^oiPi + = pn^oiPi] + ^ / mcnm-k) + y: [ $(r)^ (63) 

Z Jk „>3 n\ Jr 

The functional (158!) can be further simplified in the part of the phase diagram corre- 
sponding to the uniform phase, where (p*(r)) = Po = const, just by assuming $ = in 
Eqs.( ll02[) - (IM1) . Eq. (155i) for the uniform phase reduces to the form 

Pff\p*o) -P + PoVcoiO) + f3ff\p*o)g = 0. (64) 



D. Comparison with the Landau-type theory 

After all the assumptions and approximations described in the preceding sections, 
we finally arrived at the form of Qco, Eq. fl63l) . similar to the Landau-Ginzburg- Wilson 
and Landau-Brazovskii theories. In the original Landau-Ginzburg- Wilson and Landau- 
Brazovskii theories focusing on universal features of the order-disorder transition, it is pos- 
tulated that the effective or coarse-grained Hamiltonian in the uniform system has the form 




where 

(3H2[<I>] = \ Jj{k)C',{m~k). (66) 

The summation in (165|) is truncated at n = 6 when a tricritical point is studied. The form 
of Cl{k) is 

C',{k) = To + eok' (67) 

or 

C'M = ro + eoik-hf (68) 

in the LGW and LB theories respectively. In the original Brazovskii theory only even powers 
of the field are included {h = A3 = 0). In nonuniform systems, with the equilibrium density 
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profile p*(r), one obtains the effective Hamiltonian Heff[p*{v) + 0(r)] as a functional of the 
fluctuation 0(r). 

Note that and fl58p can be reduced to the LGW or LB functional (jHS]) with the 
following expressions for the coupling constants 

An = (3fi''\pl) (69) 

for n > 3 and 

h = f5ul!\-pi)-ii-v,o{m) (70) 

in the case of the uniform phase ($ = 0), where the explicit forms of fh^\pQ) are given in 
Appendix C for fh approximated by fl5T]) . However, in our theory (see fl59l) ) 

C:^°ik) = (3VUk) + A2 + (71) 

When Vco{k) assumes the global minimum for k = kb = 0, and can be expanded about 
k = 0, then fl^Tl) corresponds to this expansion truncated at the second order term. Trunca- 
tion of the expansion of Vco{k) is justified when the fields p*(k) with large k yield neg- 
ligible contribution to Q. This is the case when thermally excited density waves with 
large k are associated with significantly larger energy than the density waves with the 
wavenumber /c — > 0, and can be disregarded. The above conditions are satisfied when 
the global minimum is deep, and local minima, if exist, correspond to significantly larger 
values of Vco{k). More precisely, the interaction potential should satisfy the condition 
\Vcoik) - (VcoiO) + V^oK^)k^)\/\Vcoik)\ < 1 for that yield the dominant contribution to n. 
Under the above conditions our mesoscopic theory reduces to the Landau-Ginzburg- Wilson 
theory and describes phase separation. 

When Vco{k) assumes the global minimum Vco{kb) < ioi k = kf, ^ 0, and can be 
expanded about k = kb, then the truncated expansion 

VUk) = VUh) + V^o\h){k - kbf/2 + ... (72) 

yields C2°{k) similar to the LB form fIBS]) (see flTTl) . and recall that $ = in the uniform 
phase). Truncation of the expansion of Vco{k) is justified when the fields p*(k) with large 
|k — kfcl yield negligible contribution to Q. This is the case when the global minimum 
of Vco is deep, i.e. the thermally excited density waves with large |k — k^l are associated 
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with significantly larger energy than the density waves with the wavenumber k ^ kb. The 
approximate version of the mesoscopic theory reduces to the Brazovskii theory for such forms 
of the interaction potentials, and describes microsegregation, or formation of lyotropic liquid 
crystals. In the following we shall focus on systems with such forms of Vco{k) - they include 
weakly charged globular proteins, nanoparticles, colloids or rigid micells in various solvents 
that mediate effective interaction potentials. Recall that in the mesoscopic theory we choose 
the length scale R, and consequently introduce the cutoff n/R. The scale R should be such 
that the dominant contribution to Q comes from the fields p*(k) with k < ir/R, otherwise 
it should be arbitrary. The results of the Landau-Brazovskii theory describe actual ordering 
when the dominant contribution to fl depends on kb, but is independent of R. 



E. Brazovskii theory 

Let us briefly summarize the original Brazovskii theory, with the coupling constants 
expressed in terms of physical quantities according to Eqs.( [69|) - (!7T]) and (172!) . and discuss 
conditions of its validity. The Brazovskii theory is particularly simple, and analytical results 
can be obtained easily. In the original Brazovskii theory C2 is calculated to first order in 



A4, and is given by [17 1 



C2{k) = C'^\k) + ^g = T + (3%l{k - hf, (73) 



where we simplify the notation by introducing 



2Vco{kb 

and 



(74) 



T = C2ih) = c:^''ih) + ^g, (75) 

and where the dimensionless temperature is defined by 

T* = = (76) 

C2°{k) and Q are given in Eqs. (I7T]) and (l62ll respectively, and kb corresponds to the minimum 
of Vco (see the definition of Vco in Eq. (fT3!) ). Here and below k is in units and length is in 
a units. Note that (1731) is consistent with our result (|6T1) up to a correction which is of order 
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Aj. When the fluctuations with k ^ dominate, then the main contribution to Q comes 



from k ^ kb- In this case the regularized integral fl62l) can be approximated by 17|, |33| 

Q = f — ^ - / . (77) 



The integral on the RHS of fl77j) can be calculated analytically [33|. In the case of r <^ (3*V2kl, 
Q takes the asymptotic form 

Q ^r«p*.lkl + + 0(/n(7r/i?)), (78) 

where 

Q{t) = (79) 

and 

a = kl/{4nV^;). (80) 



Note the independence of the dominant term in fjTHj) on the mesoscale R. In the Brazovskii 
theory the i?-dependent terms in (ITHIl are neglected. The Brazovskii approximation 

g ^ Gir) (81) 

is valid when the second term in (ITHI) is indeed negligible compared to the first term for 
1 < R < n/kh, i.e. 

r « ktP*v;. (82) 

Note also that the original Brazovskii theory is restricted to ~ kf,, because of the approx- 
imation ( 1721) . From ( 1751) . ( IHTi) and ( 1711) we obtain the explicit expression for r = C2{kb), 

r^^ = r^^C^^h) + A,aVf* . (83) 

The grand potential functional f H71) in the Brazovskii approximation takes the explicit 
form 

/3n[p* + $] = /3n,,[p* + $] + 2aVrrv - ^^^^v, (84) 

o 

where r satisfies Eq.fl83l). and ^4 are given in (11061) and fllOSp respectively, and PQco[p* + 
$] in the Brazovskii- type approximation is given in Eg. (1631) . The second term on the RHS 



in flH^ is the explicit form of the second term on the RHS in P7|) 19l . l3l| . In calculating the 
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third term on the RHS in (H7|) . Eqs. (173!) and (l58l) . as well as the property (^(r)^""'"^)^ = 
were used. The above expression is valid provided that the condition fl82l) is satisfied. 

Global minimum of the functional fl84p with respect to $ corresponds to the stable phase 
for given pg and T*. The problem of finding the minimum of fl5^ becomes easy for periodic 
structures of given symmetry. For given symmetry $(r) can be written in the form 

$(r)=^<|.„^7„(r), (85) 

n 

where gn{'^) represent the ortho normal basis functions for the n-th shell that have a particular 
symmetry, and satisfy the normalization condition 




is the n-th amplitude. For given symmetry the problem reduces to the determination of 
the minimum of the function of variables $„• 



V. EXPLICIT RESULTS 



A. Structure of the disordered phase and the A-line 

Let us focus on the stability of the disordered phase. In the first step let us limit ourselves 
to the MF approximation and consider stability of the functional f2co- When < 0, 

then Vtco is unstable with respect to the density wave with the wavenumber k. At the 
boundary of stability the second functional derivative of VL^o vanishes for k = kb, i.e. 

Crih) = 0, (87) 

since such instability occurs at the highest temperature for given density. In the PY ap- 
proximation the above equation yields together with flTTj) and fll06p the universal curve 

where the dimensionless temperature is defined in (1761) and Vco{k) assumes the global mini- 
mum, Vco{kb) < 0, for k = kfy- Universality in this context means that the shape of Vco{k) is 
irrelevant, and the value at the minimum, Vco{kb), sets the temperature scale. For properly 
rescaled temperatures the boundaries of stability of Qco for all systems with particles having 
spherical cores coUaps onto the single master curve (|88|) . 

23 



For kb = 0, Eq. (l88i) represents the MF approximation for the spinoda' 



liquid separation, whereas for kb the above represents the A-hne 
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45l| associated with microsegregation. Similar result was obtained previously 



44j . The universal line flHHl) is shown in Fig. 3 (dashed line). 



Let us consider the actual boundary of stability of the uniform phase in the system in 
which Vco{k) assumes the global minimum for kb > 0. Beyond MF the boundary of stability 
of the grand potential Q is given by 



C2{kb) = 0. 



(89) 



In the Brazovskii theory C2{kb) = r satisfies the equation (!83|) . Eqs. (!89|) and (183!) yield 
T = independently of density. This means that when Vco{k) assumes the global minimum 
for kb > 0, the MF boundary of stability with respect to periodic ordering, Eq. (155]) . is shifted 
down to zero temperature when the mesoscopic scale fluctuations are included. At the A-line 
the structure factor S{k) = G2{k)/p*Q + 1 diverges for k ^ kb for G2{k) approximated by 
G'^{k) = 1/C2°{k). In the Brazovskii approximation the maximum of the structure factor 
is finite at the A-line, and we obtain its value from (l83l) and ([8] 



G2{kb) = 1/r 



k^Al 



(90) 



Recall that (183|) is valid provided that r satisfies the condition (!82|) . which at the A line 
takes the form 



(91) 



For Pq ^ the approximation (jHTl) is not valid at the A-line, because the LHS of ( 19111 
behaves as (po)~^ (see Appendix C), whereas the RHS is independent of pg. The RHS of 
fl9T]) depends on the system, therefore the density range for which the approximation flHTj) is 
valid at the A-line, is system-dependent. The condition (!82|) is satisfied below the A-line for 
sufficiently low T* (see fl83l) for Cco{kb) < 0). The first-order transitions in the Brazovskii 
theory occur on the low-temperature side of the A-line p/7|] , therefore the simple approximate 
theory can be applied. 

The infinite susceptibility with respect to periodic external field, resulting from singularity 
of the structure factor for finite k, is the artifact of the MF approximation. In the mesoscopic 
theory the interpretation of the A line follows from Eqs. ([9]) and (|63|1 . which show that when 
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C2°{kb) < 0, the mesoscopic density p*(r) = Pq + ^gi{r), with infinitesimal-amphtude $ and 
Qi given in Appendix D, is more probable than the mesoscopic density p*(r) = Pq. This is 
because for infinitesimal $ the remaining contributions in Eg. (1631) are irrelevant. For T* < 
T{(po) (i.e. C2°{kb) < 0) the single mesostate with periodic density, with the wavenumber 
kf, and infinitesimal amplitude, is more probable than the uniform mesostate. Recall that 
the mesoscopic density p(r) is equivalent to the set of microstates that satisfy (jl]). Thus, for 
T* < T*(po) the probability of occurrence of any microscopic state with nonuniform density 
distribution at the length scale n/khis larger than probability of finding any microscopic state 
with position-independent density (jll) on the length scale R ~ Tr/ki,. However, averaging over 
all microscopic states with different spatial inhomogeneities may lead to space independent 
average density, unless the system undergoes a first-order transition to ordered phases with 
finite amplitude $ of density oscillations. The first-order transition can be determined 
beyond the stability analysis. 

When Vco{0) < 0, then both, the A-line and the spinodal line are present and the tem- 
perature at the MF spinodal line is 

r;(p-o)=TA*(pS)^^- (92) 
^s*(Po) < T^* (pn) when Vco{k) assumes the global minimum for kf, > 0. Both lines were found 



in Ref. 



41 



42 



44| for some forms of the interaction potential. In fact there is a family of 



curves T{(pQ)\4o(^)/Ko(^fe) representing the MF instabihty with respect to density waves 
with the wavenumber k. All of them lie below the A-line when k^ > 0, or below the spinodal 
line when kb = 0. 



B. Case of weak ordering - universal features of the MF phase diagram 



Let us determine the most probable mesoscopic states for given thermodynamic vari- 
ables. The most probable mesostate corresponds to the global minimum of the functional 
PQco[Po + ^] with respect to $(r) - the position dependent density deviation from the mean 
value Pq (see ([H])). On the other hand, the minimum of Qco is equivalent to the MF approxi- 
mation, where the second contribution to the grand potential Q in Eq.f l26p is neglected. The 
minimum of the functional of the same formal structure as in Eq. (163|) was calculated in 



different context in Ref. 



13, 



19 



461 ]. and we shall not repeat the details of the calculation 
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FIG. 3: Dashed line is the universal line (j88p that represents the MF spinodal line of the gas- 
liquid separation when Vco{k) assumes the global minimum for A; = 0, or the A-line associated with 
formation of ordered periodic phases when Vco{k) assumes the global minimum for k = > 0. 
Solid line is the universal first-order transition line between the disordered and the bcc phases in 
the one-shell MF approximation. T* and are dimensionless temperature and density defined in 
Eqs. ([76]) and respectively. 

which can be found in the above papers. 

In this section we shall limit ourselves to the one-shell approximation. This approximation 
is valid when Ko(^f)) ^ Ko(^fe2), where kh2 is the wavenumber in the second shell. For 
structures possesing different symmetries the Fourier transform of gi in (l85l) has the form 



where w* is the complex conjugate of w, and ww* = 1. 2n is the number of vectors k;^ in the 
first shell of the considered structure. The forms of gi{r) in real-space representation are 
given in Appendix D for the lamellar, hexagonal, bcc and gyroid (Ia3d) structures. From 
( l63]) we obtain in the one-shell approximation 




(93) 



(3AuU^,) = /3(fi,,[p* + $i^i(r)] - Q,om)/V = lc2{k,)<l>l + + (94) 



where the geometric factors characterizing different structures are given by 




(95) 
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Note that in the one-shell approximation pAucoi^i) depends on Vco{k) only through the 
product PVco{kb) (see flTTl) ). therefore the phase diagram in variables Pq,T* is universal. 
From the extremum condition 

= C,(h)i, + + ^4? = (96) 

we obtain $i of the stable or the metastable phase. The stable phase corresponds to the 
lowest value of jSAucoi^i) for given thermodynamic variables, where $i satisfies fl96|l . The 
coexistence between different ordered phases takes place when the grand potentials for 
these phases are equal. At the coexistence of the stable ordered phase with the disordered 
(uniform) phase 

Au;,„(<l>i) = 0. (97) 

The universal MF phase diagram obtained in this way in the one-shell approximation is 
shown in Fig. 4. For ^3 7^ the ordered phase coexisting with the fluid has the bcc 
symmetry, and the uniform-bcc phase coexistence line is 



T* = (981 



as already shown by Leibler 



46|. Explicit expressions for A^ are given in Appendix C, and 
the transition line (!98ll is shown in Fig.3 together with the A-line, and in Fig. 4 together 
with transitions to the other phases. The lattice constant of the bcc phase is a = 
Along the coexistence between the uniform and the bcc phases 

^ _ JM^, (99) 

where the explicit forms of the geometric coefficients, k!'^'^ = 2/^/3 and n'l'^'^ = 15/4, were 
used. $5^^^ at the transition to the uniform phase is independent of Vco, i.e. in the MF one- 
shell approximation is universal, and is shown in Fig. 5 for our PY form of fh- For A^ = the 
transition is to the striped (lamellar) phase j46] and is continuous $1 = 0). In our PY theory 
A^ip*) = for just one density p* = p*, where p* is the maximum at the A-line, and at the 
same time the critical density of the gas-liquid separation. In the PY approximation ( ISTjl 
p* ^ 0.2457358. Note that for densities lower than the gas-liquid critical-point density p* we 
obtain a periodic array of excess number density (cluster) forming the bcc crystal, whereas 
for higher densities the bcc structure is formed by bubbles of depleted density. Note also 
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FIG. 4: MF phase diagram in the one-shell approximation. T* and Pq are dimensionless tempera- 
ture and density, Eqs. ([76|) and respectively. Solid line is the coexistence of the uniform phase 
(above in T*) and the bcc crystal. On the left from the maximum the bcc structure is formed 
by droplets (excess density), on the right by bubbles (depleted density). The bcc crystal coexists 
with the hexagonal structure along the dashed line. Again, on the left and on the right from 
the maximum the hexagonally packed cylinders consist of droplets and bubbles respectively. The 
hexagonal phase coexists with the lamellar phase along the dash-dotted line. 

that $1 ^ 1 as required by the construction of the approximate theory (see sec.IIIB). The 
difference in the average densities of the coexisting phases in this approximation vanishes. 

The diagram shown in Fig.4 is the universal 'skeleton' showing the sequence of phases: 
disordered, bcc, hexagonal, lamellar, inverted hexagonal, inverted bcc, disordered. Note the 
re-entrant melting at densities well below the close-packing density. The above sequence of 
phases agrees with the sequence observed for micellar solutions and block copolymers, and 



with the sequence of structures found in recent simulations for the SALR potential [40|, |42 |. 
However, bicontinuous cubic phases found in some of the self-assembling systems are not 
present on the universal MF diagram obtained in the one-shell approximation. 

C. Example: the SALR potential in the case of very short range of attractions 

Beyond the one-shell approximation more complex structures can be stabilized, and we 
expect 'decorations' of the universal skeleton diagram (Fig.4) with regions of stability of 
more complex structures, or with structures with large amplitudes of the density oscillations. 
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FIG. 5: The amplitude of the density profile p^{r) = + ^igi'^'^{r) in the bcc phase at the coexis- 
tence with the uniform phase (i.e., along the solid line in Fig. 4) in the one-shell MF approximation. 
$1 and the space-averaged density Pq are both dimensionless. 

However, the diagrams are no longer universal, in the sense that the stability region of more 
complex structures depends not only on the value of (3Vco{k) at the minimum aX k = kb, but 
also on (3Vco{kbn), where kbn is the wave number in the n-th shell. Therefore the details of 
the phase diagram depend on the shape of the interaction potential. Studies of the details of 
the phase diagrams in various systems go beyond the scope of this work. Just for illustration 
we quote the results obtained for the SALR potential 

W^) = -^.^2^ + A,^^H(z££) (100) 

where r and Zi are in a and 1/a units respectively, and the amplitudes are in ksTroom units. 
The SALR potential describes in particular weakly charged colloids in a presence of short- 
chain non-adsorbing polymers inducing the depletion potential, globular proteins in some 
solvents, and rigid micells. For the parameters y4„, = 140e®'^, = 30e^'^^, Zi = 8.4 and 
Z2 = 1.55, related to the colloid-polymer mixture j?!, the Fourier transform of the potential 
Vco{r) = (^{r — ^)VsALR{f) is shown in Fig.6 together with the approximation fl72l) that allows 
for the reduction to the Brazovskii theory. This potential leads to formation of small clusters, 
because kb ^ 1.94 in cr"^ units, and 2TT/kb ^ 3.24. In Fig. 7 G2"{k) = {Sco{k) - l)po, where 
Sco{k) is the structure factor in MF approximation, is shown for Po = P* ~ 0.246 and T* = 
0.15. The first peak corresponds to cluster-cluster correlations, as observed in simulations 



and experiments for the SALR systems 



39, 



42] . The chosen thermodynamic state is 
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FIG. 6: The Fourier transform Vco{k) of Vcoir) = 0{r — l)VsALRir), where VsALRir) is the SALR 



potential (fTOOj) with Aa = 140e^-^, Ar 



30e 



1.55 



8.4 and Z2 = 1.55 (solid line). Dashed line 



represents Eq. (j72p that leads to the approximate Brazovskii theory. The minimum of Vco{k) is 
assumed at /c = fcf, ~ 1.94. The second-shell value ki)2 = 2ki)/^/3 is indicated on the A;-axis. Vco{k) 
and k are in ksTroom and 1/cr units respectively. Brazovskii theory is valid when the dominant 
deviations from uniform distribution of particles correspond to k ^ k^. 



away from the A-line, and G'^{k) is a reasonable approximation for the correlation function. 
When the A-line is approached, G'^{ki)) diverges and the MF approximation fails, as discussed 
in sec.Va. In Fig. 8 G2{kb) (maximum of the structure factor) is shown for Po = p*^ 0.246 
as a function of temperature for T* < T^{p*) in the Brazovskii approximation. In this case 
the condition ( l82l) is satisfied. We verified that inclusion of the second term in Eq.f l78l) with 
R = 1 and R = n/kf, leads to G2{kf,) ^ 0.14 and 6*2(^6) ~ 0.15 respectively at the A-line, 
while G2{kb) ~ 0.16 when the second term in ( ITHl) is neglected. Note that (^2(^6) assumes 
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rather small values, in agreement with results obtained for ionic systems 

We find the diagram shown in Fig.9 by calculating (3AuJco for the gyroid phase in the 
two-shell approximation. In this phase the second shell is very close to the first shell 
{kb2 = 2ki,/\/?)), unlike in the case of the other phases (fcfe2 = V^kb, kb2 = V^h and 
kb2 = 2kb for the bcc, hexagonal and lamellar phase respectively). For this reason the sec- 
ond shell in (|85ll should be included in the case of the gyroid phase, as argued in Ref.|l9|. 
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FIG. 7: Gf (fc) = {Sco{k) - l)p^ for the SALR potential (fTOOD with Aa = UOe^-^, = 30e^-^^, 
zi = 8.4 and Z2 = 1.55 for p* = p* 0.246 and T* = 0.15. 




0.03 0.04 0.05 0.06 0.07 0.08 0.09 



FIG. 8: G2{kh) = 1/r, corresponding to the maximum of the structure factor, in the Brazovskii 
approximation ([83]) . for Pq = pi ^ 0.246 as a function of T*. 

The details of the calculation will be given elsewhere along with the results obtained within 
the present mesoscopic theory beyond MF 361]. The unit cell of the gyroid phase is shown 
in Fig. 10 for = 0.048 and T* = 0.016. Note that in this approximation we obtain for 
low enough temperatures the sequence of phases: disordered-bcc-hexagonal-gyroid-lamellar- 
inverted gyroid-inverted hexagonal -inverted bcc-disordered. Such sequence of phases is 
found in many self-assembling systems, including micellar systems and block copolymers. 
The bicontinuous cubic phase is usually found between the hexagonal and the lamellar 
phases. In simulations of the SALR potential hexagonal and lamellar phases were found 
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FIG. 9: MF phase diagram with the gyroid phase (laSd symmetry) considered in the two-sheh 
approximation, and the remaining phases in the one-shell approximation, as in Ref.[l^. The 
interaction potential has the SALR form piIDj) with Aa = 140e^-^,zi = 8.4, = 306^-^^,22 = 1-5. 
The potential is chosen for the system consisting of charged colloids in the presence of small 
polymers, similar to the system studied experimentally in Ref. [2]. The diagram is the same as in 
Fig. 5 except that the Ia3d phase is stable in the windows between the hexagonal and the lamellar 
phases. The structure of the phase in the left window is shown in Fig. 10. For more details see 



Ref. 



for similar densities as we predict 40|]. The simulations were restricted to p* < 0.25, there- 
fore the inverted structures were not found. However, for the form of the SALR potential 
corresponding to much larger Tc/ki,, spherical, cylindrical and slab-like liquid-like clusters 
were found for increasing density, and for p* > 0.35 cylindrical and next spherical bubbles 



were seen 



42| for densities that agree with our predictions. The unit cell in this system 
contains too many particles to enable observation of the ordered phases. In colloid-polymer 
mixtures spherical clusters, elongated clusters and a network forming cluster were observed 



experimentally for increasing volume fraction of colloids j?]. Relation between the experi- 
mentally observed (presumably metastable) structures and our results obtained for thermal 
equilibrium requires further studies. 

The effect of the fluctuation contribution (the two last terms in Eq. flH^ ) on the phase 



equilibria will be described in Ref. 36|], where more details on the two-shell approximation 



will also be given. In short - when the fluctuation contribution in Eq. flH^ is included, the 



32 



FIG. 10: The unit cell of the Ia3d phase in the two shell approximation for pg = 0.048 and 
T* = 0.016. In the region enclosed by the shown surface /0*(r) — Pq > 0.06. Note that two networks 
of branching regions of excess density are present in the two-shell approximation in the structure 
corresponding to the minimum of PAujco- The lattice constant of the unit cell is 27r\/6//cb ^ 8a. 



For more details see Ref. 



36|. 



stability region of the disordered phase enlarges, and the effect increases with increasing 
temperature. For low temperatures the sequence of phases is not affected by the mesoscopic 
fluctuations. 



VI. FINAL REMARKS 



We developed a formalism that has a form of the density functional theory with additional 
term in which the effect of the mesoscopic scale fluctuations is included. The additional term 
depends on the correlation function for the density-fluctuation correlations on the mesoscopic 
length scale, for which self-consistent equations are derived. The fluctuation contribution can 
be obtained when the form of the grand potential with frozen mesoscopic scale fluctuations 
is known from the microscopic theory. The theory is designed for studies of self-assembling 
systems, where nonuniform density distributions are found on the mesoscopic length scale, 
and allows for obtaining phase diagrams and structure in terms of density, temperature 
and the parameters characterizing (effective) interaction potentials, as well as microscopic 
correlation function. We present the grand-potential functional and the equations for the 
correlation functions, starting from the most general form (Eq. fl26l) ). In sec. IV we describe 
approximate version of the theory that allows for analytical calculations. The local density 
approximation was recently compared with more accurate Rosenfeld density functional for 
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the SALR potential, and good agreement was obtained for not too high densities 4^. We 
also show that our theory can be reduced to either the LGW or the LB theory. The range of 
validity of the latter approximate theories is discussed in detail. Explicit expressions for the 
coupling constants as well as the form of the grand potential and the correlation function 
are given for the Brazovskii approximation. The above mentioned equations can be applied 
to a wide class of interaction potentials describing various soft-matter systems with particles 
having spherical cores. The theory can be generalized to nonspherical cores. 

In sec.V universal properties of selfassembly are determined in the framework of the 
mesoscopic theory introduced in sec. II-IV. In addition, a particular example of the SALR 
potential is considered. We first determine the A-line and conclude that it is related to struc- 
tural changes in the disordered system. On the low-temperature side of the A-line spatial 
inhomogeneities of a particular size A appear more frequently than densities homogeneous in 
regions of linear dimension A. We obtain universal 'skeleton' diagram (Fig.4) with the uni- 
versal sequence of ordered phases that agrees with many experimental observations in a wide 
class of self-assembling systems, and with simulation studies for the SALR potential 40, l42 |. 
This universal skeleton is 'decorated' in particular systems with more complex structures, 
and modified by fluctuations, especially for high temperatures. The continuous transition to 
the lamellar phase obtained in MF becomes fluctuation-induced first order when the meso- 
scopic scale fluctuations are included. The density difference between the disordered and the 
lamellar phase and between different ordered phases vanishes on this level of approximation. 
The transitions are presumably very weakly first order, therefore it is difficult to distinguish 
them from the continuous transitions in simulations 42 1. 

By considering the gyroid phase in the two-shell approximation for a particular form of 
the SALR potential, we find that this phase is stable between the hexagonal and lamellar 
phases for low enough temperatures (Fig. 9). The bicontinuous phase can be considered as 
a regular gel. Relation between this phase and the experimentally observed gels in colloid 
-polymer mixtures j3] is an interesting question. If the gyroid phase is thermodynamically 
stable, the gel observed in experiments may result from arrested microsegregation, and its 
structure should be more regular than the structure that is formed by arrested spinodal 
decomposition into two uniform phases. In particular, triple junctions of the cylindrical 
colloidal clusters should dominate. 
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VII. APPENDICES 

A. Two-point vertex function in Fourier representation in the case of periodic 
structures 



For periodic structures with the unit cell Vu of the volume Vu, with the structure invariant 
with respect to translations by the vector P = J2i ^iPi, where are integer and the vectors 
Pi span the unit cell, we have 



C2(k,k + k')=/ / C2(ri -r2|ri)e^''-('"^-'"2)e*(^+''')""2 = (101) 

Jri Jr2 

KEe'^'^^'^'^ ''— / 6*^'^+''')"^ / e^'^-(^^-^^)C2(ri - r2|r2) = 

p Vu Jr2&Vu ^ri~r2 



5(k + k') / e^k.(r.-r.) 1 j ^^^^^ _ ^^1^^) 

■/ri-r2 Vu Jr2<^Vu 



B. Explicit forms of for periodic structures in the local density approximation 



From (|53l) we obtain 



Cr[pS, ri] = (5 fi^\pl) + /f (pS)$(ri) + 



2 ' 3! 



C2=°[p*,$,ri,r2] =/? 



Vf (pS) + /f (pS)*(ri) + CT^<^'(ri)^V(ri " ^2) + Ko(ri2) 



Cnp*,<|.,r]=/3/f (p*) + /(^)(p*)$(r 



Cr[p*,$] = /9/f (pS)- 



103) 

(104) 
(105) 



C. Expressions for the coupling constants An in the PY approximation 



The explicit forms of = fh^\po), for fh approximated by flSTj) are 



A, 



(106) 
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I2rf + 2{)rf + Sr] - 1 



(107) 



A, 



2-12r] + 30r/2 + 11277^ + A8r]'^ 

(1 - vrpf 



(108) 



D. 51 (r) for several structures 



The function gi{r) is given for the lamellar, hexagonal, bcc and gyroid phases in Eqs. fll09p - 
( 11121) respectively, with r = (xi, X2, Xa). The minimum of Ko(^) corresponds to ki,. 



g{{r) = V2cos{kbXi) 



(109) 



^?^-(r) 



cos[kbXij + 2 cos [—^ ) cos — 



(110) 



V2 



;iii) 



^Te"^ '''' 



/ kbX3 \ 



giro 



r = A^r 



'kbXi 



cos ( — ^ I 



sm 



/kbX2\ . f2kbX3 
I — — I sm ' 



V6^ 



. (2kbXi 
sm ( 



V6 



f2kbX2\ 



+ 
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